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ABSTRACT 

Context. We present a study of the Arches cluster based on XMM-Newton observations performed over the past 8 years. Unexpectedly, 

we find that the X-ray emission associated with the cluster experienced a marked brightening in March/ April 2007. 

Aims. We investigate the origin of both the X-ray continuum emission emanating from the star cluster and the flare. 

Methods. To study the time variability of the total X-ray flux, we stacked the PN and MOS data of observations performed within a 

time interval of a few days leading to the detection of the flaring episode. We then constructed two spectral datasets, one corresponding 

to the flare interval (March/ April 2007) and another to the normal quiescent state of the source. 

Results. The X-ray light curve of the Arches cluster shows, with high significance (8.6cr), a 70% increase in the X-ray emission in 
the March/April 2007 timeframe followed by a decline over the following year to the pre-flare level; the short-term duration of the 
flare is constrained to be longer than four days. The temperature and the line-of-sight column density inferred from the flare spectrum 
do not differ from those measured in the normal activity state of the cluster, suggesting that the flux enhancement is thermal in origin. 
Conclusions. We attribute the X-ray variability to in situ stellar activity: early-type stars may be responsible for the flare via wind 
collisions, whereas late-type stars may contribute by means of magnetic reconnection. These two possibilities are discussed. 



clusters is of particular interest in the context of the chemical 
enrichment of the GC region and their overall contribution to the 
high energy activity that characterises the innermost zone of the 
Galaxy. 



2. Observations and data reduction 

We selected archival XMM-Newton observations of the Galactic 
center region targeted at Sgr A*. We reprocessed data from both 
the PN and MOS cameras (IStriider et all (l200ll) . iTurner et af] 
(l200lh ) with the tasks EPPROC and EMPROC in the Science 
Analysis Software version 9.0. The main properties of the 
XMM-Newton observations employed in the present work are 
reported in Table Q] We selected good time intervals (GTI) cor- 
responding to periods of relatively low internal background. 
For this purpose, we compiled the 10-12 keV lightcurve of the 
full field of view and selected the time intervals for which the 
count rate was lower than a certain threshold. Since the selected 
observations were performed in different conditions of inter- 
nal/particle background and orbital phase, we investigated all the 
lightcurves independently and selected an appropriate threshold 
count rate for each. The thresholds used in screening the data are 
reported in TableQ] in practice, these thresholds were selected so 
as to exclude all the peaks in the 10-12 keV full-field lightcurve. 
Throughout our analysis, we have only selected single and dou- 
ble events (PATTERN <4) for the PN and up to quadruple events 
for the MOS1 and MOS2 (PATTERN<12) cameras. Similarly, 
only events tagged as real X-rays (FLAG==0) were utilised. 

It was evident that the source was seen in an anomalous state 
in the three 2007 observations (separated by four days). For this 
reason, we considered the March/ April 2007 observations (i.e., 



1. Introduction 

Three massive star clusters are located at the Galactic cen- 
ter (hereafter GC) region: the central (IRS 16) star cluster, the 
Arches cluster (hereafte r AC), and the Quintuplet cluster (e.g., 
M orris & SerabvrJ1996b. Whereas all three are ver y bright i n the 
infrar ed band (e.g. iBecklin & Neugebauerf[l975l : iNagata et al.l 
119951) . their X-ray properties differ markedly. Both IRS 16 and 
the Quintuplet cluster are not particularly bright in X-rays wit h 
L x ~10 33 erg/s (Ba ganoff et alj2003l lLaw & Yusef-Zadehl2004l) . 
In constrast, the AC has a relatively high X- ray luminosity 
(L x ~ 10 34 erg/s)(e.g. lYusef-Zadehetal.l 120021). Although the 



total mass of each of these clusters is similar at about 10 4 M 
(iFigetl 120041) . the size of the clusters are different, with the 
AC being the most compact. The X-ray diffuse emission from 
the AC has been investig ated in a number of publ i shed papers 
over t he last decade (i.e. Yusef-Zadeh et al.l 120021 : [Wang et al. 



120061: iTsuiimoto et al J 120071) . All these studies confirm the co- 
existence in the cluster of both thermal and non-thermal radia- 
tion. The former is thought to be due largely to multiple inter- 
actions between the strong winds from massive stars, with the 
high stellar concentration in the AC providin g an explanation of 
why its X-ray luminosity is enhanced (e.g. ICanto et alJ feOOO: 
Raga et al. 2001). Irrefutably, the AC, which is the densest clus- 
ter in the Galaxy with a core density of 10 s M Q /pc 3 , provides 
an excellent opportunity to study stellar processes, within the 
environment of a massive star cluster, that lead to particle ac- 
celeration and X-ray emission. The study of these massive star 
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the flare dataset) separately from the other observations. The 
region considered for the accumulation of the count rate spec- 
trum of the source is an ellipse centered at R. A. = 17:45:50. 358 
DEC=-28:49:20.08, with axes of 0.32'x0.30'. We chose an el- 
liptical region in order to enclose most of the X-ray diffuse 
emiss ion in the AC (see contours in Fig. 2 in lYusef-Zadeh et all 
120021) . Similarly a circular region of radius 0.4', centered on 
R.A.=17:45:44.081 DEC=-28:49:25.72 was used to determine 
the background. Spectra corresponding to multiple observations 
were produced by stacking the data, using the tools MATHPHA, 
ADDRMF and ADDARF. The channels in each spectrum were 
grouped with the GRPPHA tool in order to have a minimum of 
20 counts/bin, validating the use of the^ 2 statistic. 



weighted mean of the measurements (3.4+0.4xl0~ 13 erg/cm 2 /s) 
is 0.7, fully consistent with a constant background. 

We have also studied the 2-10 keV lightcurve of the indi- 
vidual March-April 2007 observations to look for evidence of 
short-term variability during the episode of enhanced activity. 
The result of this investigation is shown in Figf2] There is no 
clear evidence of a trend in this full lightcurve; a constant fit to 
the data gives a mean count rate of 0.059+0.001 cts/s, with a 
x\ rf=l-l- Although there are some low-level short-term features 
apparent in the individual observations, we conclude that these 
are most likely due to background-subtraction errors, since sim- 
ilar effects are apparent in the lightcurves of other sources in the 
field of view. We conclude that the duration of the high state 
activity is longer than ~4 days. 



3. Results 

3. 1 . The X-ray lightcurve 

To search for evidence of temporal variability in the X-ray emis- 
sion from the AC, we constructed a "long-term" X-ray lightcurve 
based on six observation epochs, which in the majority of cases 
encompassed more than one individual observation. These were 
February 2002 (1 obs), March 2004 (2 obs), August/September 
2004 (2 obs), March/April 2007 (3 obs), March 2008 (1 obs), 
and April 2009 (3 obs) - see Table Q] 

The resulting 2-10 keV lightcurve is shown in FigQ] Here 
the flux values are derived from model fits to the stacked spectra 
for each epoch using the average spectrum for the normal state 
(see below); more specifically, the points are the weighted mean 
of the PN and MOS values, corrected for absorption. A high 
state or flaring episode is evident in the 2007 measurements. As 
a check, we also considered the net source count rates measured 
in the individual cameras, which confirmed that this flaring is 
very likely a feature of the lightcurves. To quantify the statistical 
significance of this flare, we compared the weighted mean of the 
other five measurements in FigQ]of 1.53±0.07xl0~ 12 erg/cm 2 /s 
with the flux measured in March/ April 2007 of 2.4±0.2 xl0~ 12 
erg/cm 2 /s; the variability is significant at the level of 8.6cr. 
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Fig. 1. The 2-10 keV X-ray lightcurve of the AC extending 
over an eight year period. The horizontal line is the weighted 
mean excluding the March/ April 2007 flare. The red points show 
the lightcurve for the background, where the weighted mean is 
shown with a dashed red line. 



To check whether the observed flare was produced by sys- 
tematics in the background, we also built the lightcurve of the 
flux measured from the region used for the background (the 
red points in FigQ]). For these data, the x]- ec i with respect to the 
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Fig. 2. The lightcurve from the March- April 2007 observations. 
The energy range is 2-10 keV and the time binning is set to 3 ks. 



3.2. X-ray spectral analysis 

We first performed an analysis of the stacked PN and MOS X- 
ray spectra pertaining to the normal state of the AC. The spec- 
tral model comprised a collisionally ionized plasma (APEC, 
Smit h et al J 1200 lb plus a non-thermal hard tail modeled by a 
power-law continuum with a photon index T. In the present in- 
vestigation, we found it advantageous to fix T to a value of 1, 
thereby avoiding potential biases in the spectral modeling of 
multiple datasets with a different balance between the thermal 
and non-thermal contributions. This approach is consistent with 
earlie r studies that generally found T to be in the range 1-1.5 
(e.g., IWang et ail 120061: iTsuiimoto et al.1 [2007b and consistent 
with the power-law slope characterising bright 6.4 keV fluores- 
cent line-emitting regions in the vicinity of the AC (which will 
be the subject of a later study). In the analysis, w e also fixed 
the m etallicity of the thermal plasma at 2 Z Q (e.g., IWang et all 
12006b . since the spectral fitting of both the stacked and the single 
datasets did not usefully constrain this parameter. In this context, 
we note that there is strong evidence that the di ffuse plasma in 
the GC region has a super-solar metallicity (e.g., Tanaka 2002; 
iKovama et al.ll200~7h . Two gaussian lines were included to ac- 
count for the neutral (or near-neutral) Fe K a and Fe lines, 
respectively at 6.4 and 7.05 keV, the K^/Ka, flux ratio being tied 
to 0.11. Finally all the components w ere subject to interstellar 
absor ption (WABS model in XSPEC, Morrison & McCammonl 
[19831 

The results of the joint fitting of the PN and MOS spec- 
tra for the normal state of the source are reported in the sec- 
ond column of Tableland the corresponding spectra are shown 
in FigfJ] (left panel). The best fit temperature of the cluster 
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was found to be kT =1.7+0.1 keV, in good agreement with 
what previousl y measured with S uzaku dTsuiimoto et al]l2007l) 
and Chandra (Wa ng et al.l 12006). Moreover, the resulting to- 
tal 2-10 keV flux is 5.3±0.5xl0- 13 erg/cm 2 /s (1.5±0.1xl0" 12 
erg/cm 2 /s absorption corrected), translating into a luminosity of 
3.8±0.4xl0 33 erg/s. To quantify the respective contributions of 
the thermal and non-thermal emission to the spectrum, we mea- 
sured the fluxes associated with the different continuum com- 
ponents in the best fit model. For the 2-10 keV band, the split 
was 85% contributed by the thermal emission and 15% by the 
hard non-thermal tail. Both the thermal and non-thermal com- 
ponents are very likely to be manifestations of in situ activity 
within the AC. More specifically shocks produced by strong col- 
liding stellar wi nds can both heat the plasm a and supply high 
energy particles. lLaw & Yusef-Zadehl (120041) proposed that en- 
ergetic electrons within the cluster might be able to upscatter IR 
radiation from stars via Inverse Compton Effect; in this case, the 
power-law component should contribute about a sixth of the total 
flux, in excellent agreement with the measurements. 

Table 2. Results of the spectral fitting for the normal state and 
the net flare spectrum of the AC. The values are the weighted 
mean of the PN and MOS measurements. The fluxes (2-10 keV) 
are in photons/cm 2 /s. 





NORMAL 


FLARE 


N H (10 i2 cm- 2 ) 


10.4+0.9 


10.2+2.1 


kT APBC (keV) 


1.7+0.1 


1.8+0.3 


flux A/ , £C (10- 5 ) 


5.4+0.4 


5.2±0.9 




1.0 


1.0 


flu X/w (i(r 5 ) 


1.3+0.1 




fWio- 6 ) 


1.4+0.7 




flux(itr 5 ) 


7.0+0.5 


5.2+0.9 


X 2 /dof 


234/1144 


246.5/852 



Given the discovery of significant variability of the X-ray 
flux from the AC, the next step was to investigate the spectrum 
of the enhancement. For this we used the stacked 2007 dataset 
as the source event file, with the background taken to be the 
source spectrum of the AC during its normal state. The PN and 
MOS spectra corresponding to the flare (net flare spectra) are 
presented in the Figj3] (right panel), in comparison to the nor- 
mal state spectrum (left panel). The main difference between the 
spectra of the two states of the AC is in the 6.4 keV fluorescent 
line from neutral (or near-neutral) Fe; the flare spectrum does not 
show significant emission at this energy. This clearly indicates 
that the 6.4 keV Fe K a line has a different origin to the higher 
ionization 6.7 keV iron line. This observational evidence points 
to a thermal origin for the flare. This was confirmed when we 
attempted to fit the flare spectra with the same model as applied 
to the normal state data. For the flare we found that the nonther- 
mal power-law and the low-ionization state lines of Fe (at 6.4 
and 7.05 keV) were no more required. The flare spectrum is well 
described by a thermal plasma with kT=1.8±0.3 keV with the 
total flux being 3.7±0.6xl0~ 13 erg/cm 2 /s; in terms of 2-10 keV 
X-ray flux, the flare represents an increase over the normal state 
activity of about 70+13%. 

We highlight that the temperature and the N# values inferred 
for the flare spectra are in excellent agreement with the values 
derived from the normal state spectrum (see Figj3]and Table|2]). 
It would seem that the variability in the cluster is produced by an 
increase in the intensity of the thermal emission, while the other 
physical observables, namely the temperature and N# remain 



unchanged. The emission measure (EM) of this X-ray activity 
(the flare component) has been estimated from the normaliza- 
tion of the APEC component to be some 1.2x1 57 cm~ 3 . 

4. Discussion 

In this paper we have investigated the X-ray properties of the AC 
and report the first detection of significant variability in its hot 
thermal emission. The bulk of this component is thought to arise 
from the fhermalization of the strong winds expelled by massive 
stars. ICanto et alJ J2000) calculated the expected temperature of 
such a cluster wind to be 

where v M , is the terminal velocity for the stellar wind, which 
is assumed to be the same for all the stars within the cluster. 
The temperature we measured in the AC is 1 .7+0.4 keV, which 
requires a terminal velocity of 1100-1200 km/s, in good agree- 
ment with the actual wind velocities reported by ICotera et al.1 
(Il996h . In the normal state of the AC, the EM of this plasma is 
1. 4+0.3 xl0 57 ctrT 3 . If we model the AC as a sphere of 0.7 pc 
radius and set n e ~n//, the inferred electron density is 1-10 cm -3 . 
This value again compares very well to theoretical estimates ob- 
tained from modeling the X-ray emission of dense clusters o f 
massive stars (see middle panel of Fig. 3 in ICanto et al.l 12000). 
Using the superior spati al resolution of the Chandra telescope, 
lYusef-Zadeh et al] d2002l) identified three distinct spatial com- 
ponents in the AC; components Al and A2 were found to have 
strong thermal emission, whereas A3 component was found to 
be extended and more elongated towards the south-west of the 
cluster core. Overall the diffuse emission is characterized by a 
highly significant Fe K a line at 6.4 keV corresponding to (neu- 
tral or low ionization gas). Because of the lower angular resolu- 
tion of XMM-Newton, we could not resolve these three compo- 
nents and study them separate ly but simply measure a composite 
spectrum. On the same basis, iTsuiimoto et al.l J2007) have used 
Suzaku data to study a 1.4 arcmin wide region centred on the 
cluster and derive spectral parameters very similar to those re- 
ported here. 

The non-thermal radiation associated with the AC is most 
likely attributed to the presence of a sea of cosmic rays perme- 
ating the wh o le reg ion within the cluster and the surroundings; 
Valinia et al. (2000) showed that the diffuse emission from the 
Galactic Ridge can be explained in terms of Low Energy Cosmic 
Rays (LECRs) producing both the continuum and the low ion- 
ization lines from the different elements. However, the slope of 
the non-thermal continuum emission found by these authors is 
about 1.3-1.4, somewhat steeper than the value we adopt. This 
difference may reflect the extreme conditions in the vicinity of 
the AC itself. The collisions between powerful stellar winds with 
terminal velocities at or in excess of 1200 km/s can create shocks 
which accelerate CRs up to keV energies, thereby leading to a 
flattening of the spectrum in the high energy domain via either 
bremsstrahlung or inverse Compton emission. 

Here we report an unexpectedly high state of activity in the 
AC, although the XMM-Newton Off-axis PSF does not allow us 
to directly locate the flare emission wi thin the spatial substruc - 
ture observed in the AC with Chandra ( Yusef-Zad eh et al Eool . 
The spectrum of the AC in this flaring state is well modeled as a 
hot plasma with kT=l. 8+0.3 keV. This, together with the pres- 
ence in the spectrum of the high ionization lines from elements 
like S, Ar, Ca and Fe points to a stellar origin for the observed 
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Fig. 3. The 2-10 keV X-ray spectra of the AC as measured in the PN (black) and MOS (red) cameras. Left panel: The spectrum 
measured in the normal state. Right panel: The net flare spectrum measured during the flaring episode. In both cases, the best fitting 
spectral model and the corresponding residuals are also shown. 



flux increase. In the case of stellar coronal flares, the longest de- 
tected to date is nine days. For the AC we can set a lower limit 
for the duration of the enhanced luminosity of ~4 days. This 
could fit with a coronal origin for the flare. On the other hand, 
we notice that if the flare duration is a few months, its origin is 
very u nlikely to be coronal . Assuming a distance to the GC of 
8 kpc (Gilless enet al.ll2009l) . the estimated EM for the AC flare 
is 1.2xl0 57 cm' 3 , very significantly higher than values which 
are typically measured in stellar X -ray flares (~10 S5 cm' 3 ) as 



socia ted with late-type stars (e.g., Giidel 2004; iGudel & Nazel 
120091) . The X-ray flares from late type stars observed to date 
show strong correlation betwe en the peak temperature and the 
EM. lAschwanden et al.l (120081) found that the EM in a stellar 
flare is proportional to T 4 5 ; according to this relation, an EM of 
1.2xl0 57 cm -3 yields a temperature well above 10 8 K (8.6 keV), 
well beyond the value measured in the AC flare (1.8 keV corre- 
spondent to ~2xl0 7 K). We therefore dismiss the late-type star 
hypothesis for origin of the flare. 

An alternative scenario which could account for the detected 
variability entails early-type stars, possibly in binary configu- 
rations. Either an orbital modulation of the system or a sudden 
eruption by one star might produce the observed high state ac- 
tivity. In both cases the density and the velocities of the col- 
liding stellar winds might change, thus producing an enhance- 
ment in the X-ray flux. The X-ray emission arising from equally 
strong colliding winds is inversely proportional to the binary 
separation (D _1 ); for stellar parameters like the ones found in 
the AC (M ~lO~ 5 M /yr and Vco-1200 km/s) the X-ray luminos- 
ity created by colliding winds is of the order of 10 32 ' 33 erg/s, for 
D> 10-500AU. In this context, it is plausible that the enhanced X- 
ray activity might mark the p eriastron passage of a system with 
an eccentric orbit. Recently, Mau erhan et al. (120101) reported X- 
ray emission from massive stars in the GC region. Among the 
sources studied in their sample, four are located in the AC, two 
of them are inside the XMM PSF and are good candidates for 
the location of the flaring event. Another example of variability 
involving early-type stars can be found in the binary system 77 
Carinae; this has been measured to be variable with an increase 
in the X-ray luminosity of about 80%. In this supermassive col- 
liding wind binary, the X-ray luminosity reaches 10 34 erg/s and 
the flare timescales range between a few days and about 100 days 
(Moffat & Corcoranl2009l) . A likely explanation for this variabil- 
ity is the interaction of clumpy winds; this process has been ob- 
served to produced an X-ray luminosity 1.8 times the normal 



level, with a spread in the duration of the high state activity be- 
tween ~3 and 50 days, characteristics which might equally well 
describe the AC flare. 

Here we propose that the colliding wind binary hypothesis 
provides the most plausible explanation for the AC flare. On the 
basis of the lightcurve, the upper limit for the duration of the 
measured flare is about four years. Assuming a wind velocity of 
1200 km/s, the linear size of the region intersected by a propa- 
gating density enhancement over a 4-year period is about 0.005 
pc, whereas the mean distance between stars in a dense envi- 
ronment such as the AC core can be estimated to be ~ 0.1 pc. 
This argument lends further support to the stellar origin of this 
X-ray activity. Presumably, such interactions are quite probable 
in a dense cluster like the AC, especially in the southern compo- 
nent, where the highest co ncentration of massive co lliding wind 
binaries has been reported ( Yusef-Zadeh et al. 20021) . 

In the future it will be of great importance to better charac- 
terise the temporal variability exhibited by massive star clusters 
such as the AC, both from the perspective of the underlying pro- 
cesses but also in the context of the contribution of such events 
to the high energy budget of the whole GC region. 
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Table 1. The selected datasets identified by both the OBSID and the observation date. For each instrument, we report the threshold 
in the 10-12 keV lightcurve (counts/s) used for GTI selection, the total GTI exposure (ks) and the nominal observation duration (ks). 



OBSID 



Obs Date 
yyyy-mm-dd 



PN 

cut/GTI/exp 



MOS1 
cut/GTI/exp 



MOS2 
/GTI/exp 



0111350101 
0202670501 
0202670601 
0202670701 
0202670801 
0402430301 
0402430401 
0402430701 
0505670101 
0554750401 
0554750501 
0554750601 



2002-02-26 
2004-03-28 
2004-03-30 
2004-08-31 
2004-09-02 
2007-04-01 
2007-04-03 

2007- 03-30 

2008- 03-23 

2009- 04-01 
2009-04-03 
2009-04-05 



0.8/38.590/40.030 
2.0/13.320/101.170 
2.0/25.680/112.204 
1.0/59.400/127.470 
1.0/69.360/130.951 
1.5/61.465/101.319 
1.5/48.862/93.594 
1.5/32.337/32.338 
1.25/74.216/96.601 
1.0/30.114/38.034 
1.0/36.374/42.434 
1.0/28.697/32.837 



0.5/42.262/52.105 
1.0/33.070/107.784 
1.0/32.841/120.863 
0.5/80.640/132.469 
0.5/94.774/132.997 
0.8/61.002/93.947 
0.8/40.372/97.566 
0.8/26.720/33.912 
0.5/73.662/97.787 
0.5/32.567/39.614 
0.5/41.376/44.016 
0.5/37.076/38.816 



0.5/41.700/52.120 
1.0/30.049/108.572 
1.0/35.390/122.521 
0.5/84.180/132.502 
0.5/98.757/133.036 
0.8/62.987/94.022 
0.8/41.317/96.461 
0.8/27.685/33.917 
0.5/74.027/97.787 
0.5/33.802/39.619 
0.5/41.318/44.018 
0.5/36.840/38.818 



